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There have been an interesting debate on the primary source of chiral p-wave superconductivity in Sr2Ru04. 
We present a comparative study on the vortex structure between a single 2-dimensional (2D) band and quasi- ID 
band model by using Bogoliubov-de Gennes theory. The pattern of the iso-values of the local density of state 
around a vortex has a diamond shape in the quasi-lD model and is much more isotropic in the 2D model. The 
spin lattice relaxation rate well below the superconducting transition temperature is greatly enhanced in the 
vortex state in the 2D model but not in the quasi-lD model. These features can be tested by using scanning 
tunneling microscope and NMR to distinguish the models for the superconductivity in Sr2Ru04. 



I. INTRODUCTION 

The layered perovskite material Sr2Ru04 has attracted a lot 
of interests due to the experimental evidence for its spin-triplet 
superconductivity with broken time-reversal symmetry 1 . 
Shortly after the discovery of its superconductivity 2 , Rice and 
Sigrist 3 and Baskaran 4 pointed out that the superconducting 
state might be an electronic analogue of the 3 He-A phase. 
Within this scenario, assuming the simplest nearest-neighbor 
pairing interaction, the gap function in terms of the d-vector 
formalism can be expressed compactly as 5 ' 6 



d(k) = Aoz(sinfc a; + isink y 



(1) 



Sr2RuC>4 is a quasi-2D systems. Its normal state can be well 
described by multi-orbital band structure with a 2D 7 band 
derived from the Ru d xy orbital and two weakly hybridized 
quasi- ID a and j3 bands derived from Ru d xz and d yz or- 
bitals, whose Fermi surfaces are illustrated in Fig. 2(a). It 
is natural and has been generally assumed that superconduc- 
tivity arises primarily on the 2D 7 band 3 , 2D model hereafter. 
This scenario is consistent with the directional variation of the 
low temperature specific heat in a magnetic field 7 . 

However, basic questions concerning the primary source of 
the pairing remain controversial due to important discrepan- 
cies between theory and experiment, in spite of intense re- 
search work in the past almost two decades 8,9 . For example, 
the 2D model predicts an edge current in the superconduct- 
ing state, which has not been demonstrated 10 ' 11 , although this 
also has been controversial 9,1213 . Very recently, Raghu et al. 
suggested an alternative interesting possibility for the spin- 
triplet superconductivity in Sr2Ru04 12 ' 14 , where the dominant 
superconducting instability in the triplet channel occurs on the 
quasi- ID a and (3 bands, hybridized from the Ru d xz and d yz 
orbitals. This model (quasi-lD model hereafter) predicts the 
absence of spontaneous supercurrents at sample edges, con- 
trary to the single 2D band model. One expects an intrin- 
sic anomalous Hall effect due to the multi-band nature of this 
model 15 , which may explain the observation of nonzero Kerr 
effect 16 . On the other hand, the quasi- ID model would predict 
a suppression of spin density wave fluctuation at the supercon- 
ducting transition point 17 , which has not been seen in neutron 
experiments by far 18 . Therefore, it is important to explore 



the possibilities for more experimental consequences within 
the 2D band or quasi- ID band models and to test the model 
against experiments. 

In the single 2D band picture, a Cooper pair in the chiral p- 
wave state can be illustrated schematically in Fig. 1 (a), while 
its counterpart in the two quasi-lD band picture can be shown 
in Fig. 1 (b). It is highly demanding to distinguish the two 
possibilities by using available experimental probes. How- 
ever, the fragile superconducting state of Sr2Ru04 does pose 
a challenge and restriction to experimental probes. For ex- 
ample, the low transition temperature with T c w 1.5K is be- 
yond the present technical limit of using the state-of-the-art 
angle-resolved photoemission spectroscopy 19 . In order to pin 
down this controversial issue with currently available exper- 
iments, in this work we study the vortex phase of a chiral 
p-wave superconductor for both the 2D band and quasi- ID 
band scenarios. It is found that different models give rise to 
qualitatively distinguishable vortex states, a promising char- 
acteristic to identify the underlying superconducting nature of 
Sr 2 Ru0 4 . 

The rest of this paper is organized as follows. In Sec. II, 
some basic properties of the Fermi surface in the normal state 
are studied. In Sec. Ill, we construct the effective Hamilto- 
nians based on the Bogoliubov-de Gennes theory for two dif- 
ferent models as a pedagogical procedure. In Sec. IV, the re- 
sults with focus on experimental observables are presented. 
We note that the vortex state for the single-band chiral p-wave 
model has been studied in previous literature 20 . In the present 
work we use more realistic parameters for the purpose of di- 
rect comparison of two different models. Finally, a concluding 
remark is given in Sec. V. 



II. PROPERTY OF NORMAL STATE: FERMI SURFACE 

Before proceeding to the discussion of different supercon- 
ducting model, we briefly discuss the topology of the Fermi 
surface in the normal state. As we will see later, the Fermi 
surfaces of different bands give rise to different signature cor- 
responding to the superconducting state. We begin with the 
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FIG. 1. (Color online) Schematic illustrations of a Cooper pair with 
spin 5* and orbital angular momentum L for a spin-triplet supercon- 
ductor in the chiral state with d(k) formed in (a) the single band and 
(b) quasi- ID bands. In the former case, superconductivity is mainly 
from the d xy orbital, whereas the hybridization of d xz and d yz or- 



bitals play a key role in the latter case, 
distinguish these two scenarios. 



It is strongly desirable to 



Hamiltonian in the Ad — t.2 g orbital basis of the Ru ions 

g(k) o 
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where £i(fe) = — 2tcosk x — n, ^(k) = —2tcosk y — //, 
g(k) = —At' sin k x sin k y and^fc) = —2^3 (cos k x +cos k y ) — 
At' 3 cos k x cos k y — Here 1, 2, and 3 denote orbitals xz, 
yz, and xy, respectively. Hereafter we take (t,t' ,t 3 ,t' 3 ) = 
(1, 0.1, 0.8, 0.35) for the hopping parameters, and (fi, fi 3 ) are 
fine tuned such that the electron density for each band is 
equal to A/3 21 . Note that this set of parameters can repro- 
duce Fermi surface, whose shape agrees with that obtained in 
the angle-resolved photoemission spectroscopy measurement 
above T c 22 . 

After diagonalizing the Hamiltonian in Eq. 2, the Fermi sur- 
face can be obtained, as shown in Fig. 2 (a). It is noted that 
the 7 band is almost isotropic in the two-dimensional plane, 
whereas the quasi-lD a, (3 bands show significant anisotropy. 



To further illustrate this point, we show the Fermi velocity 
wV(fc) in different bands in Fig. 2 (b). We can see that the am- 
plitude of vf is maximized along the in-plane square lattice 
axis a orb for the a and f3 bands, while it is almost isotropic 
for the 7 band. 

In the subsequent sections related to the superconducting 
state, we will concentrate the discussion on the vortex state 
in order to provide a good test and the answer to the question 
in which band superconductivity takes place. As we will see 
later, the topology of Fermi surface is actually inherited in the 
superconducting state for which the band is predominant. 
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FIG. 2. (Color online) (a) Two-dimensional Fermi surface for 
Sr2Ru04 in the normal state, (b) Fermi velocity vp(k) along the 
Fermi surface line. Here Of is the angle between the point k on the 
Fermi surface and x-axis, as indicated in the figure. 



III. BOGOLIUBOV-DE GENNES THEORY OF THE 
SUPERCONDUCTING STATE 

A. Single 2D band model 

First, we consider a single 2D band model based on the d xy 
orbital. In terms of the eigen-energy E f and the quasi-particle 
amplitudes u\, v\ at the i-th site, the Bogoliubov-de Gennes 
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equations are given by 
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E 



4 - H h 



E f 



(3) 



where if. 



•j 



^3<5ij , and e is an index of the 
eigenstate. The magnetic field is introduced through the 
Peierls phase factor e lipii with tpij = A(r) • dr, where 

A(r) — y, x, 0) stands for the vector potential with mag- 
netic field H in the symmetric gauge and $o = hc/2e is the 
superconducting flux quantum. Within this choice of gauge, 
the next-nearest vortices of the square-vortex lattice are lo- 
cated at the 45° directions from the (100) direction. This 
vortex lattice configuration is suggested from the neutron- 
scattering experiment 20 ' 23 ' 24 . We have = t 3 and = t' 3 
for the nearest-neighbor and next-nearest-neighbor hopping, 
respectively. The self-consistent equation for the pairing po- 
tential is reduced to 



V 



l,j±e 



¥ u\) tanh 



2T' 



(4) 



with T the temperature. Here e = x, y, denoting the unit 
vector along the x and y direction, respectively. In this paper, 
we set V = 1.0. 

The pair potential at each site i can be decomposed into p x 
and p y components as 



A Px (r 4 ) 
Here we have denoted 



Ag i = exp 



A Sii - A_ 



(5) 
(6) 



A(r) ■ dr 



(7) 



For sinpx ± isinp^-wave superconductivity, we can define 
the paring potential as A±(n) = A Px (r») ± iA Py (r^). 



B. Quasi-ID model 

We start with an effective two-orbital Hamiltonian that 
takes into account of only the Ru d xz and d yz orbitals 12 ' 13 ' 15 . 
By assuming an effective attraction that causes the p-wave su- 
perconducting pairing, one can construct an effective model 
to study the vortex physics of the chiral p-wave superconduc- 
tors in the mixed state. Here we ignore the on-site repulsion, 
which is not expected to change our conclusions qualitatively 
in the present problem. 

After the mean-field decomposition, one arrives at the 



Bogoliubov-de Gennes equations 
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where H ijj0l p = -e^^Uj^p - 5ij5 a pn, and u 
and Vj n are the Bogoliubov quasiparticle amplitudes on the 
j-th site with corresponding eigenvalues E e . Here m and n 
denote the d xz and d yz orbitals, respectively. The hopping 
integrals are chosen as 



t a = (3 = m(n), i = j ± x(y), 

-f a^/3,i = j±{x + y), 

t' a^fi,i = j±(x-y), 

otherwise. 



(9) 



The matrix elements in the off-diagonal terms of Eq. 8 are 
obtained through the following self-consistent equations, 
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The orbital part of the order parameter can be decomposed 



A I \ Ax i 
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where we denote 

Ax i — A{ i^x,mm GXp 



Ay,i 



A_, 



(11) 
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For suvp x ± i sinpy-wave superconductivity 25 , we can de- 
fine the paring potential as 



A±(n) ee A p Jr i )±iA Py {r i ). 



(14) 



Without loss of generality, we restrict the calculations with 
A + and A_ as the major and minor components of or- 
der parameters. Two types of vortices arise depending on 
the direction of the magnetic field insofar as the chirality is 
fixed 26 . In numerical computations, the unit cell with size 
N x x N y = 33x33 and the number of such magnetic unit cells 



M x x My = 5 x 5 are used 
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C. Calculation of Experimental Observables 

In this subsection, we discuss the formalism for the lo- 
cal density of states and the nuclear spin-lattice relaxation 
rate in both models for Sr 2 Ru04. In the scanning tunneling 
microscope (STM) experiment, the tunneling conductance is 
proportional to the local density-of-states (LDOS) N(E,Ti), 
which can be calculated as 

N(E,Ti) = - ^2 l [\Ul\ 2 f'(E t -E) + \V e i \ 2 f'{E t + E)]. 

(15) 

where f'(E) is the derivative of the Fermi-Dirac distribution 
function with respect to energy. Hereafter we denote 

Ut=u\, Vt = vt (16) 

for the single-band model, and 

Ut = <ra + <,„, Vt = Km + <,n (17) 

for the quasi- ID model. 

In addition to the STM measurement, nuclear magnetic res- 
onance (NMR) is another related powerful method to identify 
distinct signatures predicted by different models. Generally 
speaking, this method is able to simultaneously shed light on 
the spatial profile of the zero-energy quasi-particles through 
the relaxation time T±. The nuclear spin-lattice relaxation rate 
we consider is given by 28,29 

i?(r,,rv) =lmx+-{r l ,r il ,iQ, n fi + ijy)/(n/T)| n _ >0 

e,e' 

X7rTf'(E e )5(E e -E e ,). (18) 

We choose r, = iv by considering that the nuclear spin-lattice 
relaxation at a local site is dominant. Then the site-dependent 
relaxation time is given by T\(r) ~1/R(r, r). Roughly speak- 
ing, T\ is proportional to the integral of the LDOS within the 
energy range < E < T. Therefore, the NMR experiment 
is also expected to provide important fingerprints for the two 
models. 



IV. RESULTS 

In this section, we will present results for the vortex states 
based on the two different models and basic methods dis- 
cussed in the previous section. 

A. Vortex Structure 

To begin with, a general picture on the vortex structure for 
the superconducting order parameters for both models will be 
shown below. We will also show the spatial dependence of 
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FIG. 3. (Color online) Vortex lattice structure in a single band model. 
Left panels are for a negative vortex lattice, and right panels for a 
positive vortex lattice. Upper panels (a) and (b): Spatial distribution 
of the major component |A+(r)|; Middle panels (c) and (d): Spatial 
distribution of the admixed component |A_(r)|; Lower panels (e) 
and (f): LDOS N(E = 0, r). The size of a magnetic unit cell is 
33 x 33. 



the pairing order parameter and LDOS at zero bias. While the 
former cannot be measured directly, the LDOS is an experi- 
mentally observable quantity in the STM measurement. 

First, we study the single band model and plot the order pa- 
rameters as a function of position in a vortex lattice in Fig. 3. 
Due to the broken time-reversal symmetry of the chiral p- 
wave state, there are two types of vortices depending on the 
direction of the magnetic field. However, the negative vor- 
tex with winding number opposite to the chirality is the stable 
state 30 . And as shown in Fig. 3, the orientation of the square 
shape is different depending on the winding. In the negative 
or positive vortex case, the shape of |A + | around the vortex 
core is nearly isotropic with minor anisotropy. However, the 
induced component |A_| in both negative and positive vor- 
tex shows similar shapes, which extends along the a-axis. 
Finally, as shown in Fig. 3(c) and (f) we also calculate the 
LDOS N(E = 0, r) which is related to the tunneling conduc- 
tance at zero bias in the STM experiment. We can see that the 
LDOS comes to a peak around the vortex core, with slightly 
anisotropic extension along the (110) direction. This is be- 
cause that the amplitude of Fermi velocity vf for the 7 band 
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is almost isotropic with slight enhancement near the (110) di- 
rection, as shown in Fig. 2 (b). Such features are qualitatively 
similar to previous results 20 . 
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FIG. 4. (Color online) Vortex lattice structure in a two band model. 
Left panels are for a negative vortex lattice, and right panels for a 
positive vortex lattice. Upper panels (a) and (b): Spatial distribution 
of the major component |A+(r)|; Middle panels (c) and (d): Spatial 
distribution of the admixed component |A-(r)|; Lower panels (e) 
and (f): LDOS at zero bias N(E = 0, r) . The size of a magnetic unit 
cell is 33 x 33. 

Now we turn to the results for the quasi- ID model. Simi- 
lar calculations as the single band case are performed with the 
corresponding results shown in Fig. 4. We find several qualita- 
tive differences resulting from a different models after careful 
comparison. First, in this model, the amplitude of the major 
component A + (r) is not substantially suppressed at the vor- 
tex core, compared with the previous model. This can be seen 
from the density scale of Fig. 4 (a) and (b) showing a variation 
from about 0.3 to 0.6, whereas the counterpart for the single 
band model in Fig. 3 (a) and (b) is from 0.1 to 0.5. Secondly, 
for the induced component A_(r), its amplitude reaches a 
maximum at the vortex core. On the contrary, |A_(r)| maxi- 
mizes near the core along the a or b axis and shows a fourfold 
symmetry. 

More importantly, one should pay attention to the LDOS at 
zero bias for this model. At a first glance, the shape of LDOS 
around the vortex resembles a rhombus ( ), with its vertices 
pointing along the a or b axis. Such a highly anisotropic 



structure is actually an indication of strong anisotropy for the 
angle-resolved Fermi velocity. This intuitive understanding is 
supported by our calculations showing that \vp \ in the normal 
state is maximized along the a or b axis for both a and j3 bands 
in Fig. 2 (b). As a result, the shape of the LDOS at zero bias, 
acting as a fingerprint of the active band where the supercon- 
ductivity arises, can be used to distinguish the two different 
models in the STM measurement. 

Moreover, we find that the peak feature at the vortex core 
is substantially smeared for the quasi- ID model (pay attention 
to the density scale of Fig. 4 (e) and (f) in comparison with the 
single-band counterpart in Fig. 3 (e) and (f)). This property is 
unique for a quasi- ID model 31 , and the underlying physics is 
the strong suppression of low-energy resonance in the vortex 
state of a quasi- ID model, which will be discussed in detail 
with more experimental signatures below. 



B. Signatures in STM and NMR Measurements 

After having a general perspective on the vortex structure 
for the two different models, below we would like to study 
and analyse the LDOS in detail. The energy dependence of 
LDOS can be directly probed by measuring the tunneling con- 
ductance with an appropriate voltage bias in the STM experi- 
ment. In Fig. 5, we show a comparison of the energy-resolved 
LDOS at various positions for the two models. In the single 
band model, we can see the LDOS at the core of the negative 
vortex has a peak at E = 0, while the peak for the positive 
vortex is slightly higher than zero. This feature, qualitatively 
consistent with previous calculations, is due to different wind- 
ing structures of the negative and positive vortices 20 . 

However, the situation is totally different in the quasi- ID 
model. Similar study on organic superconductors in a mag- 
netic field 31 indicates that vortices are strongly modified due 
to the quasi- ID nature of superconductivity. Specifically, it 
has been demonstrated that vortices in a quasi- ID supercon- 
ductor do not possess low energy mid-gap excitations. This 
extraordinary property leads to the missing of mid-gap reso- 
nance peak in the energy dependence of LDOS at the vortex 
site. To see whether this conclusion can be applied to our 
quasi-lD model for Sr2Ru04 or not, we also calculate the en- 
ergy dependence of the LDOS at different positions in Fig. 5 
(b). It can be seen that the mid-gap resonance is also absent 
in the quasi-lD model we consider here. Nevertheless, we 
want to stress that this conclusion is valid only when the inter- 
orbital hopping amplitude t' is small enough. On the contrary, 
if t' is large enough, say t' > OAt, we find that the mid-gap res- 
onance will be present again, similar to the single band model. 

Not only does the STM technique enable direct detection 
of the mid-gap excitations at the vortex core, but these ex- 
citations are also expected to be indicated in the NMR ex- 
periment. Now, we proceed to the discussion of the nuclear 
spin-lattice relaxation rate Tj~ , which can be calculated via 
Eq. 18. For the single band model, from Fig. 6 (a) we can see 
that a residual relaxation rate shows up at the vortex core as 
the temperature approaches zero, making T^ 1 deviates from 
the bulk value in the absence of a magnetic field by several or- 
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FIG. 5. (Color online) LDOS N(E, r) as a function of energy E 
at various positions for (a) the single band model and (b) quasi- ID 
model. The bulk value is obtained in the absence of a magnetic field, 
whereas "N core" ("P core") indicates the position at the negative 
(positive) vortex core. The temperature is about 0.3T C . 



FIG. 6. (Color online) The nuclear spin-lattice relaxation rate T-j - as 
a function of temperature T at various positions for (a) a single band 
model and (b) a quasi-lD model. The bulk value is obtained in the 
absence of a magnetic field, whereas "N core" ("P core") indicates 
the position at the negative (positive) vortex core. 



ders of magnitude especially at low temperatures, consistent 
with previous calculations 32 . This huge deviation stems from 
the fact that the mid gap states of the vortex core contribute 
substantially to the relaxation rate, compared to the zero field 
case. Next, we turn to the quasi-lD scenario in Fig. 6 (b). 
Note that the temperature dependence of the relaxation rate at 
different positions is almost identical in the logarithmic scale. 
The only subtle difference between the bulk and the vortex 
core is that T^ 1 is slightly enhanced at low temperatures in 
the vortex core, but the previous residual behavior is totally 
absent. This result is consistent with the previous discussions 
on the energy dependence of LDOS, which also show similar 
behaviors between the bulk and the vortex core. 

For Si"2RuC>4, the NMR measurement has been performed 
experimentally in the absence of a magnetic field down to 
0. IK 33 . And we suggest that further experiment in an external 
field should be performed as a crucial test of the two candidate 
models. 



V. SUMMARY AND CONCLUSION 

We have studied the vortex state for a chiral p-wave super- 
conductor based on a single 2D band and two quasi- ID band 
models, respectively. By comparing the two sets of results, 
we have found several distinctive characteristics in the vor- 
tex state derived from different models. Generally speaking, 
for the band structure of Sr 2 RuC>4 at the Fermi surface, the 
7 band is more or less isotropic in the x — y plane, whereas 
the quasi- ID a and j3 bands are highly anisotropic. Assuming 
that superconductivity originates from either the 7 band or a 
and j3 bands, the vortex state does inherit distinguishable fea- 
tures from the active band(s). Firstly, the shape of the LDOS 
at zero bias in the quasi- ID model shows anisotropy in accor- 
dance with the analysis on the angle-resolved Fermi velocity 
for the a and j3 bands. Secondly, the missing of the mid- 
gap resonance at the vortex core for this model leads to cor- 
responding consequences in the energy dependence of LDOS 
and the temperature dependence of nuclear spin-lattice relax- 
ation rate . All these specialities show sharp distinction 
compared with the counterparts in the single -band scenario. 
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Consequently, the STM and NMR measurements in the vor- 
tex state of Sr 2 RuC>4 are expected to be unambiguous experi- 
ments to answer the question which band the superconductiv- 
ity resides in. Before concluding, we would like to mention 
that it has been recently proposed that this disputatious issue 
related to the orbital origin of superconductivity can be settled 



by detecting the Leggett-like collective modes 
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